Monotonic growth of interlayer magnetoresistance in strong magnetic field in very 

anisotropic layered metals 
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It is shown, that the monotonic part of interlayer electronic conductivity strongly decreases in 
high magnetic field perpendicular to the conducting layers. We consider only the coherent interlayer 
tunnelling, and the obtained result strongly contradicts the standard theory. This effect appears in 
very anisotropic layered quasi-two-dimensional metals, when the interlayer transfer integral is less 
than the Landau level separation. 
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Introduction. The investigation of the angular and 
magnetic field dependence of magnetoresistance provides 
a powerful tool of studying the electronic properties of 
various metals. The Fermi surface geometry of the most 
metals has been measured using the magnetic quantum 
oscillations (MQO) of magnetoresistance. The angu- 
lar dependence of magnetoresistance also gives the im- 
portant information about the electronic structure and is 
widely used to investigate the electronic properties of lay- 
ered compounds: organic metals (see, e.g., Refs. 0-01 for 
[reviews), cuprate high-temperature superconductors. [8l-fll1| 
heterostructuresfHJ etc. 

In layered quasi-2D metals, where the interlayer transfer 
integral t z is considerably smaller than the in-plane electron 
Fermi energy, the electron dispersion is given in the tight- 
binding approximation by 



e 3D (k) « e (kx,k y ) - 2t z cos(k z d), 



(1) 



where e (k x , k y ) is the in-plane electron dispersion, k z is out- 
of-plane electron momentum, and d is the interlayer spacing. 
Xf t z still much larger than the Landau level (LL) separation 
hw c = heB/m*c, the standard theory of galvanomagnetic 
properties [l-y| works well. This theory predicts several spe- 
cial features of magnetoresistance in quasi- 2D metals: the 
angular magnetoresistance oscillations[13l Il4j and the beats 
of the amplitude of MQO. [J 

In strongly anisotropic layered quasi-2D metals, when 
~ hu c , many new qualitative effects emerge. For exam- 
ple, the slow oscillations of magnetoresistance appear [ll| [l6[ 
and the beats of MQO of transport quantities become 
shifted. [HI, [l7 1 These effects are not described by the stan- 
dard theory^ -Q because it is valid only in the lowest order 
in the parameter fno c /t z . When this parameter becomes of 
the order of unity, the standard theory is no longer applica- 
ble. 

The monotonic part of magnetoresistance also changes 
when t z < huj c . According to the standard theory, Q exter- 
nal magnetic field along the electric current leads only to 
MQO but does not influence the monotonic (background) 
part of this current. However, the monotonic increase of in- 
terlayer magnetoresistance R zz with the increase the mag- 
netic field B perpendicular to the conducting layers has been 
observed in various strongly anisotropic layered metals. 
24j This monotonic growth of magnetoresistance was at- 



tributed to the "strongly incoherent" regime, where the in- 
terlayer tunnelling described by the usual Hamiltonian term 
in Eq. ([5]) is not effective, and the new mechanisms of 
interlayer electron transport play the major role. For ex- 
ample, the variable-range electron hopping between the lo- 
calized states in strong magnetic field leads to the insulat- 
ing behavior and to the exponential dependence of inter- 
layer conductivity on temperature and magnetic field. [25j 
In another model, where the in-plane electron motion is 
fully metallic but the interlayer electron transport goes via 
rare local crystal defects (e.g., resonance impurities), the 
interlayer conductivity a zz also has metallic-type temper- 
ature dependence but decreases strongly with the increase 
of the out-of-plane component of magnetic field. 23] The 
boson-assisted interlayer tunnelling can describe only the 
unusual temperature dependence of interlayer conductivity 
at T ~ 10 — 150i\\ [26l - l28j but it does not explain its mag- 
netic field dependence. Below I show, that the monotonic 
growth of magnetoresistance R zz cx \fB~ z appears also in 
the standard model, described by the Hamiltonian in Eqs. 

in strong magnetic field at very weak interlayer cou- 
pling: fko c 3> To > t z , where Tq — H/2to is the electron 
level broadening due to impurity scattering in the absence 
of magnetic field and tq is the electron mean free time. This 
contradicts the common opinion [29j that in the "weakly in- 
coherent" regime, i.e. at To > t z , the interlayer magne- 
toresistance does not differ from the coherent almost 3D 
limit t z 3> Tq. This increase of magnetoresistance was also 
missed in Refs. 30-32], where the Born approximation has 
been incorrectly applied to describe almost the 2D electron 
system. 

The model. The electron Hamiltonian in layered com- 
pounds with small interlayer coupling contains 3 main 
terms: 



H — Hq + Hf + Hj 



(2) 



The first term Hq is the noninteracting 2D electron Hamil- 
tonian summed over all layers: 



(3) 



where {m} — {n, k y } is the set of quantum numbers of elec- 
trons in magnetic field on a 2D conducting layer, j{c m j) 
are the electron creation (annihilation) operators in the 
state {m} on the layer j, and e2D (jn) is the correspond- 
ing free electron dispersion given by 
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£ 2 d (n, k y ) = hw c (n + 1/2) . 



(4) 
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The second term in Eq. $2$ gives the coherent electron 
tunnelling between two adjacent layers: 

H t = 2t z J2 J d 2 r[*](r)*j-i(r) + *]_ 1 (r)* J -(r)], (5) 



where ^j(r) and ^(r) are the creation (annihilation) op- 
erators of an electron on the layer j at the point r. This 
interlayer tunnelling Hamiltonian is called "coherent" be- 
cause it conserves the in-layer coordinate dependence of the 
electron wave function (in other words, it conserves the in- 
plane electron momentum) after the interlayer tunnelling. 
The last term 



H, 



E 



d 3 rVi (r) ¥(r)V(r) 



(6) 



gives the electron interaction with impurity potential. The 
impurities are taken to be point-like and randomly dis- 
tributed on conducting layers with volume concentration 
rii and areal concentration Ni = riid on each layer. The 
impurity distributions on any two adjacent layers are un- 
corrected. The potential Vi (r) of any impurity located at 
point Yi is given by 



Vi (r) = US 3 (r - n) . 



(7) 



We also introduce the 2D point-like impurity potential with 
the strength Vq — U \ip (z.i)\ 2 w U/doi each impurity: 

Vi(x,y) = V S(x-Xi)6(y-yi). (8) 

In the limit, t z < To, hui c , the interlayer hopping t z must 
be considered as a perturbation for the disordered uncoupled 
stack of 2D metallic layers. The 2D metallic electron system 
in magnetic field in the point-like impurity potential has 
been extensively studied. |33l - l39| In the self-consistent single- 
site approximation the coordinate electron Green's function, 
averaged over impurity configurations, is given by 



G(ri,r 2 ,e) 



faK^n), (9) 



where VE^ fe {r\ ) are the 2D electron wave functions in 
perpendicular magnetic field, [40] and the Green's function 
G (s,n) does not depend on k y : 



G(e,n) 



fuj c (n + l/2)-S(e)' 



(10) 



where £ (e) is the electron self-energy part due to scattering 
by impurities. 

The interlayer conductivity a zz , associated with the 
Hamiltonian ([5]), can be calculated using the Kubo formula 
and the formalism, developed for the metal-insulator-metal 
junctions. [4l| In analogy to Eq. (44) of Ref. 29], 



4e 2 t 2 d 



d 2 rd 2 r' 



de 
2^ 



-n' F (s)] 



HL X Ly 

x (ImG R (r,r',j,e)lmG R (r',r,j + l,e)) 



(11) 



The angular brackets in Eq. (|11[) mean averaging over im- 
purity configurations. Since the impurity distributions on 



adjacent layers are uncorrelated, one can perform this av- 
eraging separately for each layer. The averaged Green's 
functions are translational invariant: (G R (r,r',j,e)) = 
(Gn(r — r',j,e)). Therefore, one can perform the integra- 
tion over r'j which removes the sample size L x L y : 



2a T Q 



d 2 r / de [ 



x (lmG R (r,j,e)) (lmG R (r,j + l,e)) 



(12) 



where we introduced the interlayer conductivity without 
magnetic field 



(T = e t z v 2D d/HT , 



(13) 



= 2NllI^c = m* /irh 2 is the 2D DoS at the Fermi 
level in the absence of magnetic field per two spin compo- 
nents, and Nll is the LL degeneracy per unit area. 

When the magnetic field is perpendicular to the conduct- 
ing layers, the coordinate dependence of the electron Green's 
function on the adjacent layers is the same. Then the in- 
tegration over r for the Green's function of the form ([9]) is 
very simple and gives the factor Nll: 



^Jds [-n' F (e)]^|ImGfi( e ,n)| 2 



(14) 



In the zero-temperature limit, where —n' F (e) = S(e — /J,), 
and in weak magnetic field, where the summation over n 
can be replaced by the integration over n, Eq. (|14|) gives 

*„(B)=*oTo/\lmE(ji,B)\ (15) 

in agreement with the standard theory. 

Calculation. In strong magnetic field, when haj c 3> 
T = irriiU 2 p3D — TrNiV 2 NLL/^c — fr/2i~ , one can con- 
sider each Landau level separately. In the self-consistent 
single-site approximation 33] the electron Green's function 
on each LL is given by 



c( ^ ) = g ±MzM ^E B El , (16) 

and the DoS on each LL is described by the well-known 
dome-like function 33] 



-ImG R (E,n) V \E - gQ (E 2 - E) 
= D (E) = , (17) 



2tt\E\E 9 

where the electron energy E is counted from the last occu- 
pied LL: E = e — Siu (fip, k y ) , and E g — NllVq, where the 
LL degeneracy per unit area is Nll — ^/2^l 2 Hz = eB /2ttHc. 
The boundaries of the DoS dome in Eq. ([T71) are 

E 1= E g (^--lf , E 2 = E g (^+lf , (18) 

where Ci is the dimensionless ratio of the impurity concen- 
tration to the electron concentration on one LL: 



Nt/NLL = 2irl 2 Hz md. 



(19) 



The function D (E) in Eq. (|17[) is nonzero in the interval 
< Ei < E < E 2 and normalized to unity: / D (E) dE = 1. 
The LL half-width 



T B = (E 2 - Ei) j2 = 2E g JF l cx s[B. 



(20) 
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The LL broadening Tg in Eq. (fTTj) is much larger than Tq 
and depends on magnetic field, which is emphasized by the 
subscript " B" . The ratio 



r B /r » y/4tujj c /irr > 1 



(21) 



grows as V B in high magnetic field. 

Taking zero temperature and substituting Eq. (|17[) into 
Eq. ((TU) we obtain 



er zz (£0 = 2^ 



7T 



2 LB B fl 



(22) 



where E = fi — E2,d (nF,k y ) and the real part of the square 
root must be taken, which is nonzero only in the interval 
Ei < E < E 2 - The monotonic part a zz of conductivity can 
be obtained by the averaging of Eq. (|2"2"|) over the oscillation 
period hui c : 



a zz (E) dE/?uo c 



2ctoTo 



En 



2 l^i 



1 + Ci 



In 



y/ci- 1 



When Ci 3> 1, this simplifies to 
2CTor 



7T.B, 



4T 



(23) 



(24) 



The interlayer conductivity in Eq. (j2~4")l decreases with the 
increase of magnetic field: er zz oc B~ x l 2 . Qualitatively, this 
dependence is obtained by substituting |ImE (/i, B)\ ~Yb 
and Eq. $ZT$) into Eq. ([15]): 



llmSl 



O"0 



£0 



(TO 



4fiw r 



(25) 



In Ref. [42j the qualitative arguments, similar to those in 
the derivation of Eq. (|25l) . have been applied to show the 
monotonic growth and the change in the angular depen- 
dence of interlayer magnetoresistance. However, the argu- 
ments in Eq. (l25|) are not strict, because |ImS (//, B)\ 7^ Tb, 
being a strongly oscillating function of magnetic field B and 
of Fermi level /1. Therefore, the calculated value of u zz in 
Eq. (j2~4")l is 4/-7T s» 1.27 times greater than the qualitative es- 
timate in Eq. (|25|) . and the above calculation of interlayer 
conductivity, resulting in Eq. (I24[) . is more strict than in 
Ref. H. 

Discussion. The physical origin of the decrease of the 
mean interlayer conductivity a zz can be understood as fol- 
lows. The 2D electrons in magnetic field are much stronger 
affected by the impurity potential: they become localized, 
and the energy of each localized electron state m is shifted 
by the energy W (to) ~ NiVo- This energy shift depends on 
the electron state m and on the conducting layer j. There- 
fore, when the electron tunnels between two conducting lay- 
ers, the energy of the initial and final states are different, 
which decreases the interlayer conductivity. 

The large increase of the effective imaginary part of the 
electron self energy |ImS (jx, B)\ as compared to Tq in the 



limit hjj c 3> To , i z , resulting to the decrease of the inter- 
layer conductivity according to Eq. (|2"5l) . can also be ob- 
tained by the following qualitative arguments. The average 
difference AW (to) of the energy shifts of two localized elec- 
tron states is determined by the fluctuation of the number 
of impurities effectively interacting with the localized elec- 
tron. This number is approximately Ni/Nll — c, > 1, and 
the typical fluctuation of this number is ~ \fci- The av- 
erage difference of the energy shift between two localized 



states is T 



(Awy 



w/ 



\/Tohoj c serves 



as an effective |ImE| in Eq. (|25l) . Indeed, the fluctuating 
shift of the electron energy is equivalent to the coordinate 
dependent ReS (r) in the electron Green's function in Eq. 
pTO . The averaging of the electron Green's function over 
impurity configurations is then similar to the integration 
over ReS (r) with distribution of the width Tw- For the 
Lorentzian distribution of the energy shift W 



D(W) = T w /ir {W-{W)f + T 



w 



(26) 



this immediately gives the imaginary part r^y ~ yTolI^ of 
the electron Green's function: 



(G R (e,n)) = 



dWD(W) 



E-W-s 2D (n) - iT 
1 

E-(W)- S2D (n) - i (T + T w ) 



(27) 



In the Green's function in Eq. P^|) , obtained in the self- 
consistent single-site approximation [33|, this averaging over 
the energy shifts of localized electron states is hidden, and 
the resulting value of |ImS| ~ Ty/ ^> r in the interval 
Ei < E < B2 where the DoS is nonzero. 

Eq. (f2"4")) gives the decrease of the monotonic part of con- 

_ —1/2 

ductivity a zz oc B z . It has a singularity at B = be- 
cause it is derived in the limit of strong magnetic field when 
huj c 3> Tq. In the crossover region Huj c ~ Tq > t z the above 
arguments remain qualitatively valid, but the quantitative 
dependence tr Z2 (B) requires additional calculation. 

In the calculation we assumed the normalized impurity 
concentration Cj > 1, because the numerous weak defects 
and the impurities, situated far from the conducting layers, 
are important for the lifting of LL degeneracy in all layered 
materials. [38] Therefore, Cj > 1 even in the strongest pulsed 
magnets with B - 100T. 

We do not go beyond the self-consistent single-site 
approximation [33| in studying the influence of the impurity 
potential, because further correction s g ive only the small 
tails to the DoS distribution. |34l-l36l [39l | Hence, these cor- 
rections do not change the main result. We also disregard 
the electron-electron interactions, which restricts our study 
to the limit when the Fermi energy is much greater than 
the cyclotron energy, so that many Landau levels are oc- 
cupied. The chemical potential oscillations [43l [44| are also 
neglected for two reasons: (i) they do not considerably af- 
fect the nonoscillating part of conductivity and (ii) they 
are strongly damped (almost cancelled) by the MQO of the 
sample volume. [45[ This magnetostriction was directly ob- 
served in beryllium. No chemical potential oscillations 
are observed also in very anisotropic layered organic metals 
13"- (BEDT-TTF) 2 SF 5 CH 2 CF 2 S0 3 . 



4 



To summarize, we calculate the inter layer conductiv- 
ity in strong magnetic field in very anisotropic quasi-two- 
dimensional metals. The calculation is performed in the 
framework of the coherent tunnelling model, given by the 
Hamiltonian in Eqs. ([2])-([7|). In this calculation the impu- 
rity scattering is considered in the self-consistent single-site 
approximation, which is much more accurate for layered al- 
most 2D metals than the traditionally used Born approx- 
imation. This allows to obtain the new qualitative effect: 
the strong growth of interlayer magnetoresistance with the 
increase of magnetic field along conductivity and perpendic- 



ular to the conducting layers [see Eq. (|24p]. This result may 
explain the numerous experimental observations in strongly 
anisotropic layered organic metals [D3H24T ]. where the inter- 
layer conductivity strongly decreases with the increase of 
magnetic field al ong conductivity in contrast to the stan- 
dard theory @ H ESIH- 
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